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' Abstract. In this paper, we present appendices employed in the paper "A geometric 

' level-set formulation of a plasma-sheath interface" by the authors. 
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Appendix A. Local existence of sheath solutions 



In this appendix, we present a series of a priori estimates for the approximate solutions 
constructed in Step - Step 7 in Section 7 and then give the proof Theorem 7.3. 

V 

' A.l. Basic a priori estimates. In this part, we give a priori estimates for the approximate 

'nI" . solutions constructed in Step - Step 7. 
O 

. A. 1.1. A priori estimates for Stepl . In Lemmas A1-A3, we will give a proof of the ex- 

' istence, uniqueness and regularity for n as given in Step 1 of Section 7.2. 

o ■ 

, We first consider the equation for a characteristic curve. For given (x, t), 

, In what follows, we will use calculus type estimates for the Holder seminorm. For /j G 

CO'T(A,(r; v)) i = 1,2, we have 



(A.l) 9,x(s,t,x) = v(x(s,t,x),s), x(i,t,x)=x, < s < T. 



(A.2) [[/l/2]]0,7 < [[/l]]0,7lll/2|||0 + |||/l|||0[[/2]]0,7, 

><■ (A-3) rto,7<elll^lll°[[/]]o,,. 

H ; 

_Cr ! Here [[•]]o,7 and ||| • |||o denote the Holder and esssup norms defined on the same space-time 

region. 

In the following Lemma, we use simplified notation for balls in M?: 

Bi := B{0,rb + 3Ko6*To) and B2 := B{0,rh + 6KoS*To). 

Lemma A.l. There exists a sufficiently small constant Tq > and a unique solution x to 
the equation HA . 1\) satisfying the following estimates: For < T < Tq, v G B{T), 

(1) The forward characteristic curve x(s,0,x), s > 0, x G 17^(0; v) C — JIq) hits 
the target boundary dO.^ and the ion-density in the region Ag(T; v) is given by 

n(x(t,0,x),t) =no(x)exp(-^ (V • v)(x(s, 0, x), s)ds) , xGf),HO;v). 

(2) x(s,t,x) G C^'^<{[{),T] X [0,r] X M2) and sup max \d^^x'\ < 2. 

S,t,X*J = l:2 

1 
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(3) Suppose that Vi ^ v in C^''^(A(T)) and let x.% '^'^i X the characteristic curves 
corresponding to Vj and v respectively. Then for (x, t) E Aj(T; v), 

Xi{;t,^)^x{;t,^) mCi'^([0,r]). 

(4) a(x, t) := x(0,t,x) is Lipschitz continuous in (x, t) G A(T) wi^/i a Lipschitz con- 
stant 4, i-G- 

|a(x,t) -a(y,s)| <4|(x,t)-(y,s)|. 
Proof, (i) It follows from the dissipative condition (T>2) in the definition of B(T), we have 

v(x, i) • X < - ylxp, (x, t) G (^2 - f^o) X [0, T]. 
Then we have for (x, t) G {B2 - ^q) x [0, T], 

^ 2(v(x(s,t,x),s),x(s,t,x)) < -r/o|x(s,t,x)p. 

Here (•,•) denotes the standard inner product in M^. Hence the characteristic %(s,t,x) 
satisfies 

|x(s,0,x)| < e"~|x(0,0,x)| =e~ 2 |x|. 
So x(s,0,x) has decreasing magnitude and must hit the target at some positive s. 

Let T <Tq and we define the subregions A](T; -v),Q,l{T; v) of A(r) and r2s(0) as in Step 
1 of Section 7.2.1. Then the characteristic curve x 

(s,0,x), (x, 0) G O^O) X {t = 0} hits 
the target boundary O^Iq and will provide the ion density n at the target boundary, i.e., 

n(x(t,0,x),t) = no(x)exp(-^ (V-v)(x(s,0,x),s)(is), xGJ7i(0;v). 




Figure 1 . Schematic diagram of the geometry of characteristic curves 



Remark A.l. We briefly summarize the geometry of characteristic curves in the space- 
time region A(T) (see Figure 5). 

The region Ag(T) will he completely covered by the characteristic curves x{s, 0, x), (x, 0) G 
0^(0; v) X {t = 0} and they are pointing toward the target for positive s. On the other hand, 
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all backward characteristic curves x(s,t,x),0 < s < t, (x, t) G A(r) — A^(T) will either hit 
the initial region (B{0,36*) — Oq) x {t = 0} at s = or the target boundary 90o at some 
s £ [0,t) (see Figure 5). However the latter situation will not happen, for example, suppose 
the backward characteristic curve x(s,to)Xo),0 < s < to>(xo,to) G A(T) — Aj(T) hits the 
target boundary at s = sq at yg/ 

yo := x(so,*o,xo). 

Then forward characteristic curve x{s, so,yo), s G [soj^o] have the same image of a 
trajectory as x{s,to,:x.o), s G [sq, to]i (xq, ^o) S A(T) — Aj(T), but this is impossible since by 
the strong dissipation assumption 1)2 in the Definition 7.1, no forward characteristic curves 
will be issued from the target boundary. 

(ii) The first part of the proof for (1) follows from the standard theory of ordinary differential 
equations. In fact we gain regularity in the s- variable, i.e., 

x(-,t,x)Gc2'^([o,r]). 

We differentiate ()A.1|) with respect to Xj to get 

dsd^^xHs.t,^) = Vvk{x{s,t,^),s)-d^^x{s,t,^), 0<s<T k,je {1,2}, 
dxjX^{t,t,yL) = 6jk. 

Here Sjk is a Kronecker delta function and is the k-th. component of k=l,2. 

We integrate the above equation along the characteristic curve x to see 

dxjX^'itt,^) = Sjk - Vvk{x{s,t,x),s) ■ dx^x{s,t,x)ds. 

The above relation implies 

sup max l^^ x'^l < 1 + 6Ko'^*(i - 6 sup max \dx.x''\- 

s,t,xk,j=l,2 s,t,xk,j=l,2 

Since i — ^<T^l,we have 

sup max |3s,x^l < 2. 

s,t,xk,j=l,2 

(iii) Consider the equations for Xi and X- For (x, t) G A](T; v), 

d^Xi{^,t,x.) =\ri{xM,t,^),0, f55X(^,i,x) = v(x(^,f,x),0, 

Xi{t,t,y^) = X, |x(t,t,x)=x. 

We use the above equations to calculate XiiC^ — t, x), and integrate in ^ from ^ = s 
to ^ = t to get 

Xj(s,t,x) - x(s,i,x) = (vi(Xi(^,f,x),0 - v(x(^,t,x),0)d^ 

= - r {MXii^, t, x), - Vi(x(e, t, x), O) d^ 
- f (vi(x(e, t, x), - v(x(e, t, x), O) d^. 

(A.4) 
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Here we used Xj(t,t,x) = t,x) = x and note that 
(vi(Xi(e, t, x), - Vi(x(C, t, x), O) d( 

= jj Vxv,(x(e,t,x) + si(x,(e,t,x)-x(e,i,x),e) • (x,(e,t,x)-x(e,i,x))dside 

(A.5) 

We now take the M^-norm in ()A.4|) and use HA.5|) to see 

IXj(s,*>x) - x(s,t,x)| ^ 

< 1 1 1 Vvil I |o,A(T) ^ \XiiC, t, x) - t,x)\d( + \\ \vi - v| I |o,A(T) (t-s). 

Note that Gronwah's inequaUty yields 
|Xi(s,t,x) - x(s,t,x)| 

< ll|v. - v||lo,A(T)(i - + ll|Vv,|||o,A(T)(* - s)elll^^'"W)(*-)). 

(A.6) 

By hypothesis (2) of this lemma, we have Vj — > v in C^''^(A(r)) as i — > oo, and this implies 
from (|X6|) that 

(A. 7) ||Xi(s,t,x) - x(s,t,x)||o,[o,T] ^ as i ^ cx). 

Next we show 

(A. 8) ||9sXj(-,t,x) -9^x(-,^,x)||ojo,T] ^ as i ^ oo. 

Note that ()A.1|) implies 

\\dsXi{-,t,^) - 5sX(-,*,x)||ojo,T] 

< l|vi(Xi(-,t,x),-) - Vi(x(-,t,x),-)||o,[o,r] + l|vi(x(-,t,x),-) - v(x(-,f,x), OHojo.r] 

< ll|Vvi|||o^;^(T)IIXi(->*>x) - x(-,i,x)||o,[o,T] + ll|vi - ^\\\o,A(T) ^0 as i ^ cx). 

We use ()A.8|) to show 

(A. 9) [Xi{-,t, x) - t, x)]o,^jo,T] ^0, asi^oo. 

By direct calculation we have 

\{Xi - X){si,t,x) - {Xi - X){s2,t,x)\ 

\SI-S2\^ 

< \ \dsXii-,t,x) - 9sX(-,t,x)||o,[o,T]|si - ■52^"'^ 

< \ \dsXi{-,t,x) - 9sX(-,t,x)||ojo,T]^^"'^ ^0 as i ^ cx). 

Next we show 

(A. 10) [a^Xi(-,i,x) -asX(-,t,x)]o,^,[o,T] ^ 0, asi^oo. 

It follows from HA.1|) and (|A.2|) that 
[dsXi{-,i, x) - dsX{-,t, x)]o,-yjo,T] 

< [vi(Xi(-,t,x),-) - Vi(x(-,f,x),-)]o,7,[0,T] + [Vi(x(-,t,x),-) - v(x(-,t,x),-)]o,7,[0,T] 
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< lllVv,- 



j|M0,7,A(T) 



|0,7,[0,T] 



V,' - V 



l0,7,A(T) 



as i ^ oo. 



Here we used ()A.7|) and ()A.9|) . 

Finally we combine the estimates ()A.7|) - ()A.10|) to get 

IIIXi(s,*>x) - x(s,t,x)|||i^^jo,T] ^0 as i ^ oo. 



(iv) By the triangle inequality, we have 



|a(x, t) - a(y, s)\ ^ |Q:(x,t) - a{y,t)\ ^ \oi.{y, t) - a{y, s) 



- (y,s)| |x-y| |t-s, 

Here we used (|A.7|) and hypothesis (2) of this lemma. Next observe that 
|a(x,t) -a(y,t)| = |x(0, t, x) - x(0, t, y)| 

= / agx(0,t,x + e(y-x))(iC = / VxX(0,t,x + C(y-x))-(y-x)de 
JO Jo 

< [ |VxX(0,t,x + e(y-x))||y-x|dC<2|y-x|. 
Jo 

Here we used Lemma Al (1): 

III^xX|||o,[0,T]x[0,T]xR2 < 2. 

Similarly, we have 

- Q!(y,*)|||o,[o,T]xR2 < 2|t-s|. 

Hence we have 

|«(x,t) -a(y,s)| <4|(x,t) - (y,s)|. 
Lemma A. 2. Suppose f is a scalar valued function defined on Aj(T;v) satisfying 

sup ll/(-,t)|lo,7,Qi{t;v) < OO. 

0<t<T 



□ 



Then we have 
ft 



f{x{C,0,a{^,t)),OdC 







<C,iT)l sup ||/(-,t)||o,,,oi(*;v) 



where [[ ■ ]]o7A1(T v) Holder seminorm on the space-time region, and 

Ci{T) := (t^-^ + 16^r) = 0{T^-^). 



If f is in C^''^{A.l(T;v)), then the term supQ<j<j. 1 t)| |q ,^ qi^^.^^ can be replaced by 

lo,7,Al{T;v)' ^-fi-? 

ft 



/(X(^,0,a(x,t)),0de 



0,7,Ai(T;v) 



<C^l(r)|||/||lo,7,Al(T;v)- 
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Proof. Let (x, t) and (y, s) be two points in A](T; v). Without loss of generality, we assume 
that s < t. 

Jo 0, «(x, t)),OdC - /o fixiC, 0, «(y, s)),Od( 



< 



|(x,t)-(y,s)|T 
£/(x(C,0,a(x,i)),Ode 



|(x,t)-(y,s)|7 
/o 0, a(x, t)), - 0, a(y, s)), 



(A.ll) , ,/ N / M 

^ ' |(x,t)-(y,s)|7 
The terms on the right hand side of HA.11|) can be treated as follows 

/,*/(x(e,o,«(x,t)),e)d^ 



<it-s)'-''\\\f\\kKiT-,.)<T'-' 



< ( sup [f{-,0]or(,n,iO 



■ 0<?<T 
< ( sup [f{;0] 



< 



|(x,t)-(y,s)|7 
0, a(x, t)), - 0, «(y, s)), ^ 

|x(e,0,a(x,t))-x(C,0,a(y,s))r 
2^|«(x,t)-a(y,s)r 

( sup [/(•,e)]o,,,n.(o)8^i(^'*)-(y'^)i^- 



0,Ai(T;v), 



■ 0<C<T 



Note that 

max||||/|||o Ai(^ ), ( sup [/(•, t)]o.7,f7i(t;v) ) | < sup ||/(-,t)|lo,7,Qi(t;v)- 

Hence we have the desired result. Furthermore if / is in C'^''^(Aj(T; v)), then the term 
suPo<t<T ll/(->*)llo,7,Qi(t;v) can be replaced by |||/|||o,7,ai(T;v)> i-e-, 



0,7,Ai(T;v) 



<Cl(T)|||/|||o,^,Ai(T;v)- 



□ 



Lemma A. 3. Let n he the solution of (7?) given by Then there exists a positive 

constant Ti such that n satisfies the a priori estimate: 

lll^lllo,7,Al(T;v) +max|P"n|||o,^,Ai(T;v) + IPi™lllo,7,Ai(T;v) < < T < Ti, 

where Ri is a positive constant depending on Kq,6* and 7. 

Proof, (i) Recall that n satisfies 
(A.12) 

n(x,t) = no(a(x,t))exp(-^ (V • v)(x(e, 0, a(x, t)), , for (x,t) G Al{T;v). 
Since v € ;St, we have 

ll|V-v|||o,^,A(T) < (iKo6* in A(r) 
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and hence (|A.12() implies 

I ll.nl 1 1 - <' c6TK'o(5*||„ II _ 

llFlllo,AKr;v) ^ e lFo|lo,f7i{0;v)- 

Furthermore if we assume T\ is sufficiently small enough to satisfy 
(A. 13) e6Tii^o5* < 2 

then we have 

(A-14) llhlllo,Ai(T;v) < 2||no||o,ni(0;v)- 

Next we show that no(Q!(x, t)) is in C''''''(Aj(T; v)). Let (x, t) and (y, s) be points in 
A^(T; v). Without loss of generality, we assume s <t. Then we have 



\no{cx{x,t)) - no{cx{y,s))\ 



\no{cx{x,t)) - no{cx{y,s))\ f\a{x,t) - a(y,s)|\7 



|(x,t)-(y,s)|7 |a(x,t)-a(y,s)|7 

< N]o,7,ni{0;v)4^> 



l(x,t)-(y,s) 



and hence 
(A.15) 



[N(a)]]o,7,Al(T;v) < Mo,7,ni{0;v)4^- 



On the other hand, it follows from Lemma A. 2 and the fact that v G B(T) (see {'D{2)) in 
Definition 7.1) that 



(V-v)(x(e,0,a(x,t)),Ode 



0,7,Ai(T;v)) 



< 6Ci{T)Ko6* 



We use (|X3)) and (|XT3|) to get 



(A.16) 



0,7,Ai(T;v)) 



exp - (V • v)(x(C, 0, a(x, t)), 
and then use (|X2)) . (|XT3|) . (|XT5l) and (|Xl6|) to find 

(A.17) [N]o,7,Ai(T) < 2[no]o,7,ni(0;v)4^ + 12||no||o,ni(0;v)Cim^o5*. 

Since Ci(Ti) = C'(T-|'^~'^), we have for Ti sufficiently small that 
(A.18) l2Ci{T)Ko6* < 1, T< Ti, 

so that (|A.17|) implies 

(A-19) [M]o,7,Ai(T;v) < 2[no]o,^,ni{0;v)4^ + lho|lo,ci(o)- 

Finally combine ()A.14|) and (|A.19|) to get the desired bound 

lll^lllo,7,Ai(T) ^ max{22^+\3}||no||o,^,ni(o)' < T < Ti, 
(A.20) < max{22T+i,3}5*. 



< 12Ci(T)Ko(^* 



(ii) We now need to estimate space derivatives of n. Differentiate the continuity equation 

2 

dtn + ^d^^{nvi) = 

i=l 



8 MIKHAIL FELDMAN, SEUNG-YEAL HA, AND MARSHALL SLEMROD 

with respect to Xj to find 

i=l i=l i=l 

Here ^ = dt + v -Vx- 

Integrate (|A.2H) along the characteristic curve x to obtain 

(A. 22) ~ J2 [ {'^^x.XjVi + dx.Vidx^n + dx,ndx^ViYx{LO,cx{yi,t)),^)d^. 

1=1 ''^ 

Since v satisfies 

max ^max |||5'^t'j|||o,7,Ai(T;v) + |||5"vi|||Q,j,^^(y)^ < ?>Kq5* , by the definition of B{T) 
and n satisfies 

llhlllo,Al(T;v) < 2||no||o,ni(0;v)) 

for T sufficiently small by HA.14|) . we have from ()A.22|) that 
max|||a"n|||o;^i(y ) 

< max||a"no||o,ni(0:v) +6ri^o5l|no||o.f7KO;v) + 12i^o'5*Tmax|P 

(A.23) 

We assume Ti is sufficiently small so that 

(A.24) Ko5*Ti < ^, 

then we have from ()A.24|) that 

(A.25) max|||a"n|||o,Ai(T;v) < 2max ||5"no||o,m(o;v) + ^ll"'o||om(0;v) for T < Ti. 

\a\ = l ' s\ ' / \ce\ = l J sv I / 2 > sv I y 

Next we estimate [[9"n]]o^^ ai(T;v)' = 1 using Lemma A. 2 and ()A.2|) . 

By direct calculation, we have following estimates: For T > sufficiently small, we have 
(A.26) • [[dx,no{a{x,t))]]Q^^^xi(^T.^^ < 4^^ max[5"no]o,T,,ni{0;v)> 

(A.27) . [[dl,^vM]o,,MiT;.) < 6i^o(5*)'max{22^+\3}, 

(A.28) • [[^x,v^^x^n]]Q^^^xUT■,v) < 3i^o<)'*[[Vn]]o^^^;^i{T;v) + 3i^o'^*|||Vn|||o^;^i{T;v)' 
(A.29) • [[5:r,na^,^^i]]o,7,Ai{T;v) < 3i^o<^*([[Vn]]o^^^Ai{T;v) + max|||a"n|||o^Al{T;v))• 
We combine estimates (|A.26|) - ()A.29() to get 
max[[9"n]]o,^,Ai(T;v) < 4^ max[a"no]o,7,ni{0;v) + 12^o5*Ci(r) 

\a\=l J I J s V 1 / \a\=l 

(A.30) X (5* max{227+i, 3} + niax[[5-n]]o,^,Ai(T;v) + max |p-n|||o,Ai(T;v) 
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We assume Ti is sufficiently small so that 

(A.31) i^o<J*Ci(ri) < 1 

and hence for T G (0, Ti], HA.30|) implies 

max[[a"n]]o,^,AKT;v) < 2^7+1 max[9"no]o,^,ni(0;v) + max{24^ 3} 

(A-32) + 2max||9"no||o,ni(0;v) + ^ll'^o|lo,nl(0;v)• 

We combine ()A.25|) and ()A.32|) to obtain 

(A.33) max|||9"n|||o^^^Ai(^.^) < 25* max{22T+i, 4}. 

(iii) Now we estimate the time derivative of n. Recall that n satisfies 
(A.34) dtn + V ■ (nv) = 0. 

Next we use (|A3f)j) . (jAl33|) and (|X34|) to see 

ll^*"'llo,7,Ai(T;v) ^ 2 max ||c?"n||o^^^jYi(r;v) ■ lll'^lllo,7,A(r) + ll"'llo,7,Ai(T;v))ll^ " '^llo,7,A(r)) 
(A.35) < 15Ko(<5*)^max{22T+i,4}. 

Finally we set 

Ri{KQ,5\-f) := max{22T+i,3}5* + 25*max{22T+i,4} + 15i<Co(5*)^ max{22T+i, 4} 
to see that ()A.20p . (|A.33|) and (|A.35|) imply the desired result. □ 

A. 1.2. A priori estimates for Step 2. In this part, we will give existence, uniqueness and 
regularity for the function ( as given in Step 2 of Section 7.2. 

Recall from Step 2, C satisfies the exterior Neumann problem for Laplace's equation at 
given time t £ [0,T]: 

AC(-,t) = 0, xGJli, 



(A.36) 

Lemma A. 4. ho € C^''^{di}o x [0,T]) and satisfies 



■ 1^0 = ho, X G dQo and lim VC = 0. 

|a;|^oo 



lll^o|||i,7,anox[o,T] < (J* + 6CoKoRi{S*)'^ , 
where Co is a positive constant. 
Proof. Recall that ho = dtg — (nv) • Uq- Then 

lll^o|||l,7,9Cox[0,T] < \\\dt9\\\ 
(A. 37) + |||rzi?iZ^oi|||l,7,Of^ox[0,T] + |||™2^^02|||l,7,anox[0,T]- 
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Since the product of Holder continuous functions is again Holder continuous (see [31], pg- 
53), we have 

llh^l^^Ol|||l,7,9nox[0,T] + lll"'^2^'02|||l,7,aQox[0,T] 

< Co 1 1 1 n| 1 1 X [0,T] 1 1 1 ^1 1 1 1 1,-7, dQo X [0,T] 1 1 1 1^01 1 1 1 1,7,9^0 x [0,T] 

+ 1 1 1 "-1 1 1 l,7,9r2o X [0,T] 1 1 1 ^2 1 1 1 l,7,9r2o x [0,T] 1 1 1 ^^02 1 1 1 l,7,9Qo x [0,T] ) 

(A.38) <6CoKQRi{6*f. 

Here Cq is a positive constant and we used Lemma A. 3, max | ||t;j| ||]^ ^ ^^^^-j < 3Kq5* , and 
inequalities 

lll^llll,7,9nox[0,T] < lll"-||ll,7,Ai(T;v) ^nd 1 1 1 1 1 l,7,9f^o x [0,T] < 1 1 l^^il I ll,7,A(T) • 

Hence in (|A.37|) . we use (|A.38|) and the assumption (A2) of Section 7.2 to get 
(A.39) |||/io|||i,7,anox[o,T] < S* + 6CoKoRiiS* f. 

□ 

Recall the annulus region (??) in Section 7.2: 

= {jceR^ : ^ <\^\ < 26*}. 

The following existence and uniqueness result of two-dimensional exterior Neumann prob- 
lem (|A.36() is due to the results of Bers and Finn-Gilbarg |3Uj . 

Lemma A. 5. Suppose the boundary data ho is in C^''^{di^Q x [0, oo)) as provided by Lemma 
A. 4. Then there exists a unique solution C up to constant of IIA.S6\} satisfying the following 
estimates: For the compactly supported subset fi^. ofQi, we have 

(1) 

ll|VCIIIi,7,Q,x[o,T]+ sup max|P'^C(-,t)lllo,7,n. <^i> < t < T, 

0<t<T \a\=3 

where Ri is a positive constant which depends only on 6* . 
(2) Let /iq"^ G C^''^{d^o x [0, 00)) be a sequence of boundary data satisfying the bound 
1[TM) and 

/iq"^ ho in C^''^{dQo x [0, 00)) as n ^ 00. 

Suppose VC^"^ and VC are the corresponding solutions to the above exterior Neu- 
mann problem \A . 30^) for data /iq"^ , h^ respectively. Then we have 

^C^ri) _^ (C7i'^(n* X [0, T])f as 00. 

Proof, (i) Since is compactly supported in l^i, it follows from the interior Schauder 
estimates {^^, Section 6.1), we have 

l|VC(-,i)lli,7,n,+maxP"C(-,t)llo,7,n. 

a|=3 

< Co(rJo, f^*)|||/io|||l,7,aOox[0,T] 
(A.40) <Co{no,n^){6* + 6CoKoRi{6*)^), 0<t<T. 
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Let < s < t. Then it follows from Laplace's equation and the boundary condition that 
A( ^("' l t'g"-^^ ) =0, xeJ^i, 0<s<t<T, 

We now apply the global Schauder estimate ([SI], Section 6.2) to get 
|VC(x,t)-VC(x,g)| 

\t - s\-y 

< Ci(r2o, r2*)|||/io|||i,7,anox[o,T] 

(A.41) <Ci{no,n^){6* +6CoKoRi{S*f), x G f^i. 

We take the supremum over t 7^ s to get 

(A.42) sup I^C(x,^) - VC(x,g)| ^ Ciino,n,)i5* + 6CoKoRii6*f), x G fii. 

Let (x, t) ^ (y, s) and without loss of generality, assume that x / y, t / s. From ()A.40|) 
and ()A.43|) . the Holder quotient satisfies 

|VC(x,t)-VC(y,g)| ^ |VC(x,t)-VC(x,g)| |VC(x,g)-VC(y,5)| 
|(x,t)-(y,s)|7 - \t-s\^ |x-y|^ 

< [VC(x, ■)]o,^,[o,T] + [VC(-, s)]o,7,Q4s_) 

< (Co(0o,0,) + Ci(J7o,a))(5* + 6Coi^o^i('^*)'). 
Taking sup over the time-space region $7* x [0, T], we have 

(A.43) [VC]o,7,n.x[o,T] < (Coi^o, + CiiQo, ^*)W + 6CoKoRi{5*f). 

Similarly we can estimate max |||9"Clllo,7,c, x[o,T] to g^t 

(A.44) max|||5"C||lo,7,fl.x[o,T] < C2{no,n,){6* + 6CoKoRi{6*f). 

\a\—2 

Finally we combine ()A.40|) . (|A.4.S|) and (|A.44|) to get 

ll|VCIIIi,7,n,x[o,T] + sup max|||a"C(-,i)lllo,7,n* < ^1, 

0<t<T |a!|=3 

where Ri := (Co(17o, + Cii^o, Q^) + Cii^o, n^)){d* + 6CoKoRi{d*f). 

(ii) The difference ^^"^ — C satisfies 
('A(C(")(-,t)-C(-,t))=0, xGJ^i, 

I V(C^"^ - C) • i^o = 4"^ - /lo, X G aJ^o and C^"^ - C = on (9B(0, 35*). 
By the Schauder estimates (Section 6.2 in |S1]), we have 

\\V&\-,t)-VC{-Mi,r,n,<C\\hP-ho\h,r,dno- 
Letting n — > cxd, it follows from the above inequality and hypothesis (2) of this lemma that 
(A.45) VC^")(-,i) ^ VC(-,t) in (C^'^(f]i))2, 0<t<T. 

For the time-estimates we apply the same method as in (i) to get 
(A.46) VC(")(x,.)^VC(x,-) in(C7l'-([0,T]))^ x G J].. 
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We combine (|A.45|) and (|A.46|) to see 

V^W ^ VC in {C^'^ni X [0, T]))2. 

This yields 

VC(")^VC in X [0,T]))2. 

□ 

A. 1.3. A priori estimates for Step 3. We now need to give existence, uniqueness and 
regularity for the interface of Step 3 of Section 7.2. 

Lemma A. 6. 1. Assume that (" satisfies estimate (1) of Lemma A. 5: 

l|VCIIi,7,n.x[o,T] + sup maxP°C(-,i)|lo,7,n. < ^1- 

0<t<T|a|=3 

Then there exists a unique solution for the interface system (11 ) satisfying 

{e,ns,r) G (C7i'^(M x [0,T]))3 and {dpe,dpns,dpr){-,t) G t G [0,T]. 

Moreover, we have 

• ll^lll,7,Rx[0,T] + ll"'s||l,7,Mx[0,T] + I kl ll,7,Mx[0,T] 

< 2^||0o||l,7,R + 1 1 "-50111,7,18 + lko||l,7,R) i 

• ||9/36'(-,t)||2,7,R + ||<9/3ns(-,t)||2,7,M + P/jr (• , t) 1 1 2,7,R 

< 2(^||5/36'o||2,7,M + P/3nso||2,7,R + P/jT'ol |2,7,m) • 

2. Let 

V((") ^ VC in (C^'^(ll, X [0,T]))2 as given by (2) of Lemma A.4, 

and let (0*^"^ r and {9,ns,r) be the solutions of the sheath system corresponding to 
VC^"^ and VC respectively. Then we have 

(0("\n("\r(")) ^ (0,n„r) m {C^'^R x [0,T])f as n ^ oo. 

Proof. The result is just continuity with respect to data for the hyperbolic system (??). 
The proof of convergence in C^(M x [0, T]) follows from the argument in |25j . The proof of 
Holder norms (^^'^(M x [0,r]) is similar to that of (see gOl. □ 

A. 1.4. A priori estimates for Step 4- We next present the existence, uniqueness and 
regularity of the ion density n in the region A^(T) given by Step 4 of Section 7.2. 

Lemma A. 7. Let n be the ion-density obtained from Step 4- Then the formulas in Step 4, 
namely n satisfies the differential equations: 

' d 

—x{s,to,a)=v{xis,to,a),s), s > to, 
a 

— \nn{x{s,tQ,a),s) = -{V ■ w){x{s,tQ,cx),s), 
K as 



subject to initial and boundary data: 

x{to,to,cy.) = cx and n{a,tQ) 



no{a.) to = 0, 
ns{a,to) to > 0, 
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are indeed valid. Furthermore for sufficiently small T , the following estimates hold: 

n G C"^''^(A^(r; v)) and 1 |i,'y;A2{T;v) < -^2, 
where R2 is a positive constant depending only on Kq,S*,'J. 

Proof. The proof follows from Remark A.l, i.e. since backward characteristics starting at 
a point (x, E Ag(r; v) can be traced back to a point (q:,0) in the absence of the sheath 
interface, the presence of the sheath interface means backward characteristics must hit either 
a point in r2^(0; v) or a point in the sheath interface. Furthermore, the segment of backwards 
characteristic between (x, t) and {a, 0) can hit the sheath interface at most once. Indeed, the 
backwards characteristic can enter but not exit the domain As(T; v) = Al{T;v) n A^(T; v) 
through the interface surface at time < t < T. This is because initially uq = —v on 5(0), 
hence |v • + 1| < e on S{t) for < t < T by (Al) and Theorem ??, where T > is 
sufficiently small depending only on initial data, boundary data, and e > 0. By choosing 

e > sufficiently small, the vector field for the characteristic — = v(x,f) always points 

dt 

into the domain As{T; v) at the point of intersection with the sheath interface. 

Hence the formulas follow from (??). Furthermore the regularity estimates in the state- 
ment of the lemma can be obtained in a similar manner as in Lemma A. 3. □ 

We combine Lemma A. 3 and Lemma A. 7 to get the regularity result for n in the sheath 
region. 

Lemma A. 8. For sufficiently small T, we have 

neC^'^{As{T;v)) and ||n||i^^^A4r;v) < «i + ^2- 

A. 1.5. A priori estimates for Step 5. We next give the existence, uniqueness and regu- 
larity for the function (j) defined in Step 5. 

Consider Poisson's equation on the space-time sheath region As(T;v): Let t £ [0,T] be 
given and (f> satisfy 



(A.47) 



A0 = n in v), 

V</> • = g on Sr^o a-iid (p = —InUs on S{t). 



Lemma A. 9. Let n he an ion density in the sheath region As{T; v) and satisfy the a priori 
estimate in Lemma A. 8. Then Poisson's equation \A.4'T^ has a unique solution cj) satisfying 
the following estimate: 

ma,x |||9°(/>|||o,^,A,(r;v) + sup max (/.(•, t))||o f^ (4 ) < i?3. 

l<|o|<2 o<t<r|al=3 

Here is a positive constant only depending on KQ,5^i,i = 1,2 and 5* respectively. 

Proof, (i) Differentiation of ()A.47|) with respect to t shows that dt4> satisfies the mixed 
Dirichlet-Neumann problem for Poisson's equation. 

= — div(nv) inr2s(i;v), 

1^0 = dtg on 5i7o and dtcp = —dtlnris on 
where we used V0 ■ v = and Vn • = on the interface S{t), <t <T. 
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By the direct application of Holder estimates of the first derivatives given in ( j34j , Section 
8), we have 

Here Ci depends on and S{t), but we can choose uniform Ci independent of t and 
depending only on 5if and 5* for sufficiently small T, < t < T. 
On the other hand, since v S B{T), we have 

• I l^t^l li 7 ns(t-v) — ^* assumption (A2) in Section 7.2, 

where C2 and C3 are some positive constants. 
It follows from the interface equation (??) that 

dtUs / 4 sin j3 sin 96 \ / 2 sin /3y cos ' 



/4smpsmfc'y\„ ^ /2smpl/ cost/\ „ 
\ r J \ rn, / 



Us \ r / \ rUg 

^_ ^ VC • (cos6',sin6') 

We use the above relation and the estimates from Lemma A. 5 (1) to obtain 

||(9tlnns||2,7,M < C{6^i,6*). 
Here C{5^i,6*) is a positive constant depending only on 5^i,6*. Hence we have 
(A.48) snp \\dtcl)i;t)\\^^u.^.< R^^oiKo, 6,1,6,2,6*) forte [0,r]. 

0<i<T 

(ii) It follows from the Schauder estimates (Section 6.2 in [31]) that 

ll'/'(-'*)ll2,7,n4t;v) < C'4(||n(-,t)||o,^,Q4i;v) + I Isl" , 1 1 1,7,9^0 + 1 1 In n<, (•, t) 1 12,7,1 

(A.49) < R3,i{Ko,6,i,6,2,6*), 0<t<T. 

Here C4 depends only on the and S{t), but again we can choose C4 depending only on 
6, and 6* for sufficiently small T,0 < t < T. 

Let (x, t) and (y, s) be any points in As(T;v). Without loss of generality, we assume 
that < s < t. By assumption (A4) of Section 7.2, we have a contracting interface so that 

flsit; v) C flsis; v), 0<s<t<T<^l. 

Hence x G v). Then inequality HA.48|) implies, for x G fls{t;v) 

(A.50) max | |a"(/<(x, •)! |o,7,[o,T] < ^3,o(i^o, '^^i, <5*2, 5*)t1-^. 

l<|o|<2 ^ ' 

We combine ()A.49|) and ()A.50|) and choose T sufficiently small to get 
(A. 51) ^max J||a°(/>|||o,^,A4T;v) < R3,2{Ko,6*i,6,2,6*). 

(In fact the above argument holds for the expanding interfaces as well). 

(iii) On the other hand, dxi4>,i = 1,2 satisfies 
\ Adx^cf) = dxiU inl7s(t;v), 

1 Vdx,(t) ■ vq = dxiQ on d^lo and dxiCf) = -d^^ Inug on S{t). 



(A.52) 
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Again, it follows from the Poisson equation and the Schauder estimates (Section 6.2 in |34p 
that 

Px, 1 1 2,7,n4i;v) 

(A.53) < C5(p,,n(-,i)||o,^,o.(t;v) + \\dccM-,mi,^,dno + Px. lnn,(-,t)||2,^,, 

Here C5 depends on S{t), but we can choose depending only on = 1,2 and 6* for 
sufficiently small T,0<t<T. 

The first two terms in the right hand side of HA.53|) can be bounded by a quantity 
depending on 5* using Lemma A. 7 and assumptions (A1)-(A2) of the boundary data in 
Section 7.2, i.e., 

(A.54) l|9x,n(-,t)||o,^,n4t;v) + t)l|i,7,9no < C^. 

Here Cq is a positive constant depending only on and 6* . 

Now we estimate the third term H^j.. lnns||2,-).,R as follows. It follows from (??) that we 
have 

2 sin /3 

Oxi = 5/3 

r 

and similarly we can express 8x2 in terms of 9/3. Therefore we have 



(A.55) ||9,,lnn,(-,t)|||2, 



dxins{-,t) 



ns{-,t) 



<Cj i = l,2. 

2,7,R 



Here C7 is a positive constant depending only on S^i and 5* . 
Combining estimates HA.54() and HA.55|) . we obtain 

maxP"0(-,t)||2,^,f^^(i.,) < R3,3{Ko, 6,1, 6,2,6*) for t G [0,r]. 

|a|=l "' ^ ' ' 

The above inequality implies 

sup max||a"0(-,t)||o,^,o.(t;v) < R3,3{Ko,6.i,6,2,6*) for t G [0,r]. 

In particular we have 

(A.56) sup max||aX-,t)||o <i?3,3(^o,'^*i,'^*2,'5*) forte [0,r]. 

We set RsiKo, 6,1, 6,26*) := R3,2{Ko,6,i,6,2,6*) + R3,3{^o,S*i,^*2,6*) and use (|A.53|) and 
HA.56|) to get the desired result. □ 

A. 1.6. A priori estimates for Step 6. In this part, we give the existence, uniqueness and 
regularity for the ion velocity u defined in Step 6 of Section 7.2. 

Consider the Burgers' equation with a known source V0: 
(A.57) ^tu + (u • V)u = V0, (x,t)G]R2^M+. 

Lemma A. 10. Suppose the source V0 satisfies the estimates obtained in Lemma A. 9. Also 
assume initial data uq satisfy the assumption (A3) of Section 7.2 so that 

(1) V(t> G Ci'^(A,(r; v)) and V0(-, t) G C^^^i^sit, v)); 

(2) for each ct G M^, the real parts of the eigenvalues o/Vuo(q!) are non-negative; 
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the there is a positive constant T2 such that \A.57^ has a unique solution u G C'^''^{As{T;v)) 
satisfying 

detr(Q:,t) > and 

(A.58) u{xit,0,a)) = uo{a)+ [ V4>{x{s,0,a),s)ds, t e [0,T2], 

Jo 

u(x,t) •x< -^|x|2, (x,t) G iB{0,n + 6Ko6*T2)-no) x [^Ts], 

where T2 is a positive constant and 

^^^^'J^'^^ = u{x{t, 0, cx),t) and r(a, t) = Vu(a, t). 

Proof, (i) Along the particle path 0,0;), system ()A.57() becomes 

(A.59) — = Vd), where — = 5* + u • V 

^ ' Dt Dt 

Any smooth solution of ()A.58|) will satisfy 

(A.60) =^Hx{t,0,a),ty, x(a,0)=a, ' ' = uo(a). 

Since V(/>(-,i) is Lipschitz continuous and uniformly bounded, there exists a unique charac- 
teristic curve 0, a) satisfying (|A.60() locally in time t. Now we integrate ()A.60|) to get 

(A.61) dx{t, 0,cx) ^ ^ /■* 

dt Jo 
and integration of the above equation yields 

X{t,0,a) = x(0,0,q) + tuo(a) + / / V(^(x(s, 0, a), s)dsdti 

Jo Jo 

(A. 62) = a + tuo{a)+ / V(^(x(s, 0, ct), s)(is(iti. 

Jo 

Next we differentiate ()A.62|) with respect to a to get 

(A. 63) T{a,t) = I + tVuo{(x) + [ [ {V (E)V)(t){x{s,0,a),s)T{a,s)dsdti. 

Jo Jo 



Set 

y{t) = \r{a,t)-I-tVuo{a)\, 
where | • | denotes any norm on 2 x 2 matrices so that we have 



Since 



y{t)< [ [ ' \iV ^V)(l){xis,0,a),s)\\T{a,s)\dsdti. 
Jo Jo 



\T{a,t)\ = |r(Q;,t) -/-tVuo(Q:)| + |/ + tVuo(Q;) 
< yit) + l+dit, 
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where sup |Vuo(q;)| < di, we have from (|A.63() that 

y{t)< f f ' d2{y{s) + 1 + dis)dsdti, 
Jo Jo 

where |V V(/)(x(s, 0, a)| < d2 for a € ri(0), < s < T. Hence we have 

y{t) <d2[^+ dijj + d2 y{s)dsdti 

and by Appendix B 

y{t) < dsf on < t < T, 

for sufficiently small T, i.e., 

|r(Q;,t) -I-tVuo(Q;)| < dgt^ 0<t<T. 

Define 

tD{t, a) := r{a, t) - I - tVuo(Q;). 

Then we have 

\D{t,OL)\ < d^t for some constant da > 

and 

det(r(Q;, t)) = det(/ + t(Vuo(Q;) + D{t, a))). 

Let Ai(Q:, t) and Xi{a),i = 1, 2 be the eigenvalues of a matrix Vuo(q!) + D{t, a) 
and Vuo(q:) respectively. Then we can see that 

Xi{a,t) = Xi{(x) + 0{t). 

By the Cayley-Hamilton theorem, we have 

det(r(Q, t)) = {tXi{a, t) + l){tX2{cx, t) + 1) 

= {tXi{a) + l){tX2{a) + l) + 0{t^) 
= det{I + tVuo{cx)) + 0{t^). 

As long as < t < T ^ 1, the sign of det(r(Q;, t)) will be determined by det ^/ + tVuo(Q:)^ . 

Next we calculate det(/ + tVuo(Q:)). Let us set the characteristic polynomial of Vuo(q:) 
by P(q, A). Then we have 

P(a, A) = det(Vuo(Q) - XI) = (Ai(a) - A)(A2(q) - A), 

where Aj(a) are the eigenvalues of Vuo(q;). Hence 

det (^I + tVuo(Q:)) = f det (Vuo(a) + t~^I^ = t'^P{a, 

= t^iXi{a) + t~^){X2{a)+t-^) 
= {tXi{a) + l){tX2{a) + l). 

Since by assumption (2) above, real parts of the eigenvalues of the Jacobian matrix Vuo(f>:) 
are nonnegative, we have 

det r(Q, t) > n^=i [1 + t Re Xg{a)] + ©(t^) > 0, < t < T <Cl. 

Hence the Lagrangian map is a C^-diffeomorphism locally in time. 
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(ii) It follows from (|A.58|) that 

u{x{t, 0, a),t) = uo{a) + / V<f>{x{s, 0, a),s)ds. 

Jo 

We take an inner product with 0, a),t) to get 
{u{x{t,0,a),t),x{t,0,ot)) 

= (uo(a), x{t, 0, a)) + / {V^{x{s, 0, a), x{t, 0, a))ds 

Jo 

(A.64) = (uo(a), x{t, 0, a) - a) + (uo(a), a) + / (V0(x(s, 0, a), 0, 



Since 



\x{t,0,a) - a\ = j v(x(s,0,a),s)(ii 
and |||V(A|||o,A4T;!) <«3, 



Hence in (|A.64I) . we have 

(u(x(a,i),t),x(a,t)) < (!>Ko{5*fT - |a| |2 + iJgT. 
Now we choose T sufficiently small so that 

Then we have 



On the other hand, since 



^|x(0,0,a)|2 < _2r,o|x(0,0,a)|2 = -27?o|ap, 

we have 

0,0)1 < |q;| for t < T < 1. 

Therefore we obtain 

(u(x(t,0,a),t),x(t,0,a)) < 0, a)|2. 



□ 



A. 1.7. A priori estimates for Step 7 of Section 7.2. Finally we prove the existence of 
a linear extension map and some estimates of the extension. 

Lemma A. 11. Let S{t),t E [0, T] be the C'^''^ -regular simple closed convex curve in M? 
provided by Lemma A. 6 such that S{t) lies inside the annulus fi^, and v) is the corre- 
sponding sheath region < t < T , T sufficiently small. Then there exists a bounded linear 
operator £{■ -,5 (t)) : C2'T(0,(t; v)) ^ C^'^^fti) satisfying 

(a) £{u\S{t)) = u, inils(i;v), 

(b) £{u\S{t)) has support in B{0,35*), 

(c) |||^^(u|5(t))|||2,^,nj < ^o|||u|||2,^,o,(t;v), 

where Kq is independent of t £ [0,T] and il* is the annulus region 

Proof. Since the proof is rather long, we delay its proof until Appendix C. □ 
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A. 2. Continuity of T . In this part, we establish the continuity of T which in turn imply 
the existence of a fixed point of J- . 

Lemma A. 12. Let j he a continuous function such that 

f{t)<Co + Ci{f{t))\ t>0, 

/(O) < Co and CoCi < ^, 

o 

where Cq and Ci are positive constants independent oft. Then we have 

fit) < 2Co. 

Proof. Define 

F{k) = Cik^ -k + Co. 
Then by direct calculation, we have 

-I + 2C0C1 , 1 

ToinF^k) = — < at A; 



2Ci 2Ci 

Now we denote ri and r2 by the roots of F{k) = such that ri < r2. Then by direct 
calculation, the smallest root ri satisfies 

On the other hand, since F{f{t)) > 0, we have two cases: 

either / < ri or / > ^2, 
however since /(O) < Cq <ri and f{t) is continuous, we have 

/ < ri < 2Co. 

□ 

Proposition A.l. There exists a positive constant T such that the map J- with .F(v) := 
£{xi;S{t)) is a well-defined map from B{T) to B{T). 

Proof. For the time being, we assume T sufficiently small so that 
(A.65) T < min{Ti,T2}. 

So all estimates in the previous lemmas hold. 

(i) By the construction of u in the sheath region As(T;v), we have from solving (|A.57() 
along the characteristic 

fA66) l'uo(a(x,t)) + /o V0(x(s,t,x),s)ds to = 0, 

^' ^ \-i^iaix,t)) + JlV(l)ixis,t,K),s)ds to>0. 

In (|A.66|) . we have 

\Wio\\o,n4o) + i|||V(/'|||o,A,(r;v) *o = 0, 

l + (t-to)|||V0|||o,A4T;v) to>0. 

This of course implies 

ll^i|lo,A,(T;v) < lkio|lo,n,(0) + ^lll^'/'lllo,A,(r;v) 

(A. 67) < ||^ioi|lo,n,(o) + by Lemma A.8. 



\'^i\\o,AsiT;v) ^ 
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On the other hand, we use Lemma A. 2 to obtain 

(A.68) [[«i]]o,7,A4T;v) < ^4,7,^40) +Ci(r)i?3. 

We combine ()A.67|) and (|A.68|) to get 

(A.69) l|willo,7,A4T;v) < ll^oi|lo,7,n.(o) + TRs + Ci{T)Rs. 

We assume T sufficiently smah so that 

(A. 70) TR; + CiiT)R;<j. 

Here we used Ci(r) = OiT^-"'). Hence we have from (TOol) that 

6* 

(A.71) max||ui|||o^^^;^^(^.^) < max | |no»| |o,^,n,(o) + y- 

(ii) We differentiate the momentum equation 

dtU + u-Vu = V(j) 

with respect to to find 



(A.72) D{^^ ^ j^^^^^^^^^^^^^^^ ^ d,,{d,M 

where = dt + u -V. We integrate (|A.72|) along the characteristic to get 



(i) if to = 0, 

dc,^Ui{x,t) = dx^uoi{a{x,t)) -'^ / (^{dx^Uj){d^^Ujj^{x{s,t,x),s)ds 
Jo 

(ii) and if Iq > 0, 

dx^Uiix,t) = -dx,^i^i{a{x,t)) / (('9a;fc'Uj)(5x,'Uj))(x(s,i,x),s)ds 



(A.73) + / d.^^{d^^^){xis,t,x),s)ds. 

Jto 

The above equalities yield 

max max 1 1 13" Mil I In A (t-v^ < maxmax ||9"iiodln o m-i + 2T( max max || |9" {14 1110 a (t-v 

+ rmax|||9°</)|||o_A.(T;v), < t < T. 

Since T <C 1, it follows from Lemma A. 12 that 

maxmax|||a"Mi|||o,A4r;v) < 2( max max || |5"uoi||o,Q4o) + ^max 1 1|9"(/>|| Iq^A^t^^^ 

z=i,2\a\=l \t=L,z\a\=l |c*l=2 ^ ' '/ 

(A.74) < 2maxmaxP"noi||o,j^,(o)+2ri?3. 

«=i,2 a =1 
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On the other hand, it fohows from the inequahties (|A.73() and (|A.2() that 



maxmax[[9°'Ui]]o,^,A4T;v) 

< maxmax[(9°'uoi]o^o (0) + 4Ci(r) (2 max max | jS'^-uoil |o 0^(0) + 2TR^ 

4=1,2 |a|=l V j=l,2 |al=l 

X maxmax[[9"ni]]o,^,A^(T;v) + Ci{T)R^. 



i=l,2 |al=l 

We assume that 



4Ci(r)(2maxmax||a"txoi||o,A.(T;v) + 2ri?3) < \- 

\ 1=1,2, |a|=l / Z 



Here we used Ci{T) = 0{T^-^). 
Then we have 

(A.75) maxmax[[a°u,]]o,^,A4T:v) < 2 max niax[[Vuoi]]o.^,n,(o) + 2Ci{T)R^. 

We combine ()A.74|) and ()A.75|) to get 

maxmax |||9"iij|||o^ A (T-v) < 2maxmax ||9"uoj||o7 n m + '^TR^ + 2Ci{T)R^ 

i=l,2 ja|=l ' '' ^ ' 1=1,2 |o|=l ' 

We assume again that T is sufficiently small so that 

(A.76) 2ri?3 + 2Ci(r)i?3 < y- 

Then we have 

5* 

(A.77) maxlllV-Uilllo^^^A^Tiv) < 2 max HV-uoil 1 0,^,^40) + y- 

(iii) We differentiate ()A.72|) with respect to xi to obtain 

(A.78) =dl^^{d,,^). 
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We integrate the equation HA.78|) along the characteristic curve to find 
(i) if to = 0, 



j=i Jo 

+ [ 5x'fc5^,(4','A)(x(s,i,x),s)ds; 
Jo 

(a) if to > 0, 

d^^^dr„Ui{x,t) = -dl^^^Ui{cx{x,t)) 
+ I 5xfeX,(9x.'A)(x(s,i,x),s)ds. 

J to 

(A.79) 

Then it fohows from HA.79|) that 

maxmax|||5°Ui|||o,A4T;v) 

1=1,2 |«|=2 

/ 5* 
< niaxmax||5"%||o,n,(o) + 6T(2niaxmax||5"noi||o,-y,^ + - 

1=1,2 |a|=2 V 4=1,2 |o|=l 6 

X maxmax |||9"u, |||n A i't.^,\ A-TR-x. 

i=l,2|a|=2 «NI0,A4T,v) -i 

We choose T sufficiently small so that 

5*\ I 6* 
!maxmaxP"noi||o,T,,f24o) + yj < 2 TR3 < —. 

Then we have 



5* 

(A.81) maxmax|P"ni|||o,A4T;v) < 2maxmax ||a°noi||o,o4o) + 

We need to check the Holder seminorm of d"ui. Again we use (|A.79|1 to find 

max max[[9°ni]]o,^,A4T;v) 

1=1,2 \a\=2 

< rnax[[5"Mod]o,7,A,(T;v) + 6C'i(T)(2maxmax||5"noi||o,^^^^ + — 

/ 5* 
xmaxmax[[a"ni]]o A4T;v) + 6C'i(r)l2maxm^^ + — 

1=1,2 |o|=2 \ i=l,2 |a| = l o 



Here we have used (IA.2D. 



X (2inaxmax ||9"uoi|lo,n4o) + If) + C'i(^)^3 
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We assume that T sufficiently is sufficiently small that 
(A.82) 6Ci(T)(2inaxmaxP°noi||o,^,n.(o) + j) ^ ^ 

6Ci(T)(2maxmax||5"noi||o,-y,n4o) + 
(A.83) x(2maxmax||9"%||o,n4o) + ^) +C'i(r)i?3 < ^. 

Hence we have 

5* 

(A.84) maxmax[[5"ui]]o ) < 2maxmax[a°uo4 n.(o) + Tf- 

i=l,2\a\=2 i=l,2\a\=2 

We combine (|A.81|) and (|A.84|) to get 

S* 

(A.85) maxmax|||5"Mi|||o,^,A^(T;v) < 2maxmax P°Mo»||o,7,n,(o) + 

1 — i-')'^ |ci| — 2 1- — -1-1-^ \(y.\ — 2 O 

We combine ()A.71|) . ()A.77|) and ()A.85|) to get 

(A.86) max V max IP^-Ujlllo a.(T;v) < 2max V max ||9"Moi|lo,7,n4o) + < 3(5*- 

«=i,^ — \a\=k i=i,2 — \a\=k 

(iv) It follows from the Burgers' equation ()A.57|) that 

max \\\dtUi\\ Io,7,a.(T;v) < 2 max \\\ui\\ |o,7,a4T;v)) ( max max 1 1 I^^Ui 1 1 |o,^,a,(T;v)) 

+ ll|V</'|||o,7,A.(T;v) 

< I8{6*f + R3. 

(A.87) 

Finally we combine all estimates (|A.86|) and HA.87|) to get 



max ( max IP^Milllo^ Am ) < 35* and max|||9(Mi|||Q Am < ( 18(5*)^ + -Rs )• 

i=l,2 \ |a|<2 / 2=1,2 wi V y \ 

By the construction of extension of u, we find 



(A. 88) (a) max ( max |||(?"uj|||Q^(y) I < 3Ko6* , 

i=l,2 \ |a|<2 / 

(A.89) (b) max|||9iMi|||o,^^A(T) < Ko(l8{6*)^ + R3 

Here we notice that the norm || • ||o,7,ns(t;v) m Appendix C can be generalized to the space- 
time norm ||| • |||o,7,a(T)- 

Finally the estimates l|X88)) and (|X89)) show that u G B{T). □ 
We set 

As{T;r) : the sheath region determined by the interface r , 
and recall that an interface S{t) is represented by the radial function r{-,t). 
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Lemma A.13. ([S]) Let r < 7, 

C^'^{R X [o,r]) as given by Lemma A. 5 and 
Uj G C^''^(As(T; Tj)) : be associated solutions of |^ . 5? j /or each i as given by Lemma A. 10. 

Let £{ui{-,t);ri{t)) be the extension o/uj(-,t) with 

f(u,(-,t);ri(t)) tnC''^{A{T)). 

Then we have 

w = £(w ). 

As{T;r) 

Proof. The proof follows from a straightforward modification of the proof in as hence 
is omitted. □ 

Proof of Theorem 7.3 

Let {vj} be a convergent sequence in Bt in the topology of T (see (??)) such that 

Vi ^ V in C^'^(A(r)) and ^"vi ^ ^"v, in C°'^ \a\ = 2, < r < 7. 

By Proposition A.l, ^(vj) is well-defined as an element of B{T) for each i and the sequence 
{^(vj)} is uniformly bounded in T. Since the Arzela-Ascoli theorem implies the compact 
imbedding oi C^'^f {A{T)) into C^'''{A{T)) with < r < 7, we have a convergent subsequence 
which we still denote by (vj,jr(vj)): 

^(vi) ^ w in C7^'-(A(r)). 

We claim: 

(A.90) J"(v) = w. 

Proof of the claim: Let (Xj, n^, Ui,Si, (pi) and n, u, S, (p) be the quantities corresponding 
to Vj and v respectively. 

Step L Suppose that 

Vi ^ V in Ci'^(A(r)) as i ^ 00. 
Then it follows from Lemma A.l (2) that Xj(-,i,x) x(-,t,x) in C^''''([0, T]) as i ^ 00 
and hence since v G C^''^(A(T)), we have 

V-Vi(x,(e,t,x),e) ^ V-v(x(C,t,x),0 in Ci'^(af]o X [0,r]) asi^oo. 

We use Lemma A. 2 to get 
(A.91) 

/ V-Vi(x,(C,t,x),Ode^ / V-v(x(^,t,x),e)(ie inC7i'^(aJ^oX [0,r]) asi^oo. 
Jo Jo 

On the other hand, since Qj — > ct in C^''^(9$7o x [0,T]) as i — > 00, we have 

(A. 92) no^cti (x, f)) ^ no(a(x, t)) in C7i'^(ar?o x [0, T]) as i ^ 00. 

Here we used the fact that uq is in C^''^(M?). Recall the formula for n: 

n(x, t) = no(a(x, t)) exp ( - ^ (V • v)(x(C, 0, q(x, t)), 0^^?^ 

We now combine (|A.9H) and HA.92|) and the above formula to see 

nj(x, t) ^ n(x, t) in C^''^(9ilo X [0, T]) asi^oo, 



APPENDIX TO PAPER; A GEOMETRIC LEVEL-SET FORMULATION OF A PLASMA-SHEATH INTERFACE 

which in turn imphes 

hoi = dtg - (njVj) ■ v ^ ho = dtg - (nv) • v, in C^'^((9r2o x [0, T]). 
By Lemma A. 5 and Lemma A. 6, we have 

(A.93) VQ^VC in C^'^(j7, X [0,r]) as i ^ cx), 

(A.94) {0^,ri,ni) ^ (0, r, n) in C^'^(R x [0, T]). 

^tep //. Let As(T;v) be the sheath region determined by v. Since ^(v) is uniquely 
determined by v on the sheath region, once we can show w satisfies equations (??)-(??) 
and the interface conditions: 

(A.95) u = -iy and V(/> • = on 

for the orthogonal flow in the sheath region As(T; v), we will have 

w = ^(v) inAs(T;v). 

So let us proceed to show that w satisfies the sheath system (??) and boundary conditions 
HA.95|) in As(T;v). Let O be any open set compactly supported in As(T;v). Then by 
TKm^ and (f04|l . since ^ r in C^'^{R x [0,r]), we have 

OcAsiT;vi) i>N. 

For i> N,we know that {ni,Vi, (j)i,Ui) satisfy 

A(f)i = m, 

dtUi + (uj • V)ui = V(j)i, 
and 

{ni,Vi,4>i,Ui) (n,v, (/),w) in C^'^(d5), 
and hence we find in the limit as i ^ oo, 

'atn + V- (nv) = 0, (x,i)GC', 
A(j) = n, 

dtw + (w • V)w = ^(f). 

Next we check the boundary conditions on the sheath interface. Since by (|A.94|) 6i ^ 9 in 
C^''^{M.) and i^i = (cos sin we obtain 

I/, ^ ly, in C^''^(R x [0,r]). 

On the other hand, we have 

Vcpi ■ Vi = and Uj = —v>i on 5j. 

Letting i — > oo, we see 

V<j) ■ v = and w = —u on S. 

Hence we have shown that w satisfies the sheath system (??) in the sheath region and 
boundary conditions HA.95|) . By the uniqueness of the construction, we have 

J^{v) = w onAs(T;v). 

Step III. Recall by (lOIl) that we have 

ri^r in C^'^(M X [0,T]) and J^(vj) ^ w. 
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Then by Lemma A. 12, we have 

w = Sivf ) = JF(v). 

A,(r;v)^ 

Hence we showed that J- is continuous in the C^''^-topology. Since ^ is a continuous map 
on the compact and convex set B{T) of C^'^ space, by the Schauder fixed point theorem, 
has a fixed point u such that 

T{u) = u. 

This u is a desired smooth solution of the sheath system. This completes the proof. 
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Appendix B. Gronwall-Bellman type inequality 

In this appendix, we prove the Gronwall-Bellman type inequality. 

Let / be a real valued positive continuous function and suppose a nonnegative real valued 
function y satisfies the following integral inequality: 

y{t)<fit) + c^ f r y{T)dTdh. 
Jo Jo 

Then y satisfies 

y(.t) < f{t) + c^ [ [' f{s)exp[c{t-2ti + s)]dsdh 
Jo Jo 

= f{t) + ( max f{T)]o{t^), ast^O. 

Vre[o,t] / 

Proof. Let us set 

w(.t)= I r y{T)dTdti. 
Jo Jo 

Then we have 

(B.96) y{t)<f{t)+^w{t). 
It is easy to see that 

w"{t)=y{t), u;(0)=0, w'{0)=0. 
In (|(].97|) . we have a differential inequality for w: 

w"{t) < c^w{t) + f{t). 
Now we introduce another dependent variable u defined by 

w{t) = exp(ct)ti(t). 
By direct calculation, we obtain a differential inequality for u: 

u" + 2cu' < /(t)exp(-ct). 
We multiply an integrating factor exp(2ci) to get 

exp(2ct)'u')' < f{t)exp{ct). 



Next we integrate the above inequality to get 

ft fti 

u{t) < 



[ [' f{s)exp[-c{2ti-s)]dsdti, 
Jo Jo 



where we used 

u{0) = 0, u'(0) = 0. 

This implies 

wit) = exp{ct)u{t) 

Ct fti 



< [ [ f{s) exp[c{t - 2ti + s)]dsdti. 
Jo Jo 
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In (Bl), we have 

y{t) < f{t) + f [ f{s) exp[c(t - 2ti + s)]dsdti 

Jo Jo ^ ^ 

< fit) + c^( max /(r)) / / exp[c(t - 2ti + s)]dsdti 

\T&[0,t] ) Jo Jo 



f(t)+( max /(t)) Oif) ast^O, 

Vr6[0,t] / 



where we used 

ft 



r f 1 r 1 

/ / exp[c(t-2ti + s)]dsdii = ^ - 1 + -(e-'^* + e^*] 
Jo Jo C I 2 



as 
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Appendix C. Extension Theorem 

In this part, we present an extension theorem for C^''*'-functions defined on the sheath 
region Qs{t]^),t £ [0,T] to the bigger domain Qi := 5(0,35*) - Qq. 

We first consider an upper bounds for the length of a convex polygon and a simple closed 
convex curve inside the annulus A{ri,r2) defined by 

A{ri,r2) := {x G : ri < |x| < r2}. 

Lemma C.l. Let V and C be a convex n-polygon and a convex curve inside the annulus 
A.{ri,r2) respectively. Then we have 

1{V) < 2iTr2 and 1{C) < 27rr2, 

where 1{V) and 1{C) denote the lengths of the polygon V and the curve C respectively. 

Proof, (i) Let V = 'P(xi, • • • , x„) be a convex n-polygon whose vertices are xi, • • • , x„. We 
choose any point cq inside V, and we set 

yj : the intersection point with a ray cqXi and a circle B{0,r2). 

Then it is easy to see that 

(C.97) l{V{^i,--- ,x„)) </(P(yi,--- ,y„)). 

On the other hand we know that 

(C.98) /(7^(yi,--- ,yn))</(5(0,r2)) = 27rr2. 

We combine ()C.97|) and (|C.98|) to obtain 

1{V) < 27rr2. 

(ii) Let C be a simple closed convex curve lying inside A(ri,r2). Note that for any simple 
closed convex curve there exists some polygon whose sides are parts of supporting lines of 
the given convex curve. Choose a sufficiently small positive constant tq > 0. Since the 
curve C is compact, there exists a finite open cover of C consisting of balls with a center Xj 
and a radius ro, say, 

C C U*^ii?(xj, ro), where Xj G C. 

Consider an M-polygon consisting of parts of supporting lines at Xj,i = ,M and 

denote it by "P. Then it follows from the result of (i) that 

1{C) <1{V)< 2^r2. 

□ 

Next we present the existence of a continuous linear extension operator from C^'^{ris{t; v)) 
to C^'^{VLi). Even though the construction of this extension operator can be found in the 
literature, see for example H, 28, 34, we slightly modify the proofs given in books P1 0HII51] 
for our purpose. 

The proof of Lemma A. 10: We first consider the local extension near one generic 
point on the interface and then glue these local extensions together using the standard 
partition of unity to get a global extension. Let t be given. 
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Step I (local extension): Let xq be any generic point on the interface S{t). Then there are 
two cases: either S{t) is flat near xq, lying in the plane or it is not flat near xq. 

Case 1: S{t) is flat near xq lying on some line. 

For simplicity, we assume xq = (01,02) and the plane is {x2 = 02}- We choose an open 
ball i?(xo,'r) such that 

{b+ := 5(xo,r) n {X2 > 02} C B{0,36*) - !^,(t; v), 
I B- := 5(xo,r) n {xa < 02} C ns{t;v). 



Let / be any C^''*'-function defined on Qsii'i v). We extend / to the ball B^ U B as follows. 

f{xi,X2) :-- 



6/(xi, 2a2 - X2) - 8/(xi, 802 - 2x2) + 3/(xi, 4a2 - 8x2), if (xi, X2) G B~^, 
/(xi,X2), if (xi,X2) e B~. 



Notice this choice of / is not the same as given by Evans ^Hl) since he only desired C^- 
regularity. We have used a special case of the result given in JT1 . 

We claim: / is C^'-^ in the bah B. 

We need to show all partial derivatives are continuous at xq = (01,02). Let us write 

/~ := / , := f\B+- By direct calculation we obtain 

• d^J~{xi,X2) 

= Gd^Jixi, 2a2 - X2) - 8d^Jixi,3a2 - 2x2) + 3a^J(xi, 4a2 - 3x2), A; = 0, 1, 2, 

• dx^f~{xi,X2) 

= -Qdx2f{xi,2a2 - X2) + 169^.2 /(xi,3a2 - 2x2) - 95x2 / (2^1, 4a2 - 3^2), 

• dlj~ {X1,X2) 

= Qdlj{xu 2a2 - X2) - 32a2^/(xi, 3a2 - 2x2) + 275^ J(xi, 4a2 - 8x2), 

• dx^dx:,f'{xi,X2) 

= -6dx^dx2f{xi,2a2 - X2) + I6dxidx2fixi,3a2 - 2x2) - ^d^^dxj(xx,\a2 - 8x2) 

Now evaluate the above identities on the line {x2 = 02} to see that extended function / is 
in the ball B and we have 



and [C/~]o,7,B+ ^ ^^'^\Qlj\^n,B-- 
Hence we have 

ll/ll2,7,B<151||/||i,^,S-. 

Case 2: S{t) is not flat near xq. 

Since the interface S{t) is C^''''-regular, we can find a C^''''-mapping $ with inverse 
such that <I> straightens out S{t) near xq. We write y = <I>(x),/'(y) := /(<l>~^(y)). We 
choose a small ball B as before. Then as in Case 1, we extend /' from B^ to B and get 

llfll2,7,B<151||/'||2,^S-. 
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Let W := <^-^{B) and := Then we have 

ll/lk7,W'<151||/||2,-y,Ty-- 

Now we glue local extensions together using the partition of unity to get a global extension. 

Step II (Global extension): We will extend / defined on Q,s{t;w) to the bigger domain 
such that the extended / has support in Ui. Let ri be a sufficiently small number satisfying 

< ri < min((5*,0.5(5* 



Then for such ri, we choose points Xj(i = 1, • • • ,M{t)) on the curve S{t) such that neigh- 
boring Xj's are located by the part of curve with length r except one pair of points, i.e., 

/(part of an interface curve connecting Xj and Xj+i) = ro, i = 1, • • • , M{t) — 1, 
Z(part of an interface curve connecting yiM{t) ^^'^ xi) < r. 

Then the number M{t) of such points are bounded by 

-liSit))- 



M{t) < 



n 



+ 1, 



where the bracket is the greatest integer function. Then by Lemma C.l, we know that 



M{t) < 



n 



+ 1, te[0,T,]. 



As in Step I, we extend / to -B(xj, r) for each i, and denote /j by the extended function. Now 
take an open set Wq whose closure is a compact subset of fis(t) and Os(t) C Wo{t) U ^U^q^ 

Wi{t)^ . Let {ki} be a partition of unity corresponding to the open covering {Wi(i)}^Q-' of 
fls{t;v) and define 

M(t) 



i=0 



It follows from Step / that 



We take Kq to be 151 



ri 



+ 2 ) to obtain the desired result. 
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